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We present an intermediate coupling scenario together with a model analytic solution where the 
non-Fermi-liquid behavior in the underdoped cuprates emerges through the mechanism of Fermi 
surface (FS) reconstruction. Even though the fluctuation spectrum remains nearly isotropic, FS 
reconstruction driven by a density wave order breaks the lattice symmetry and induces a strong 
momentum dependence in the self-energy. As the doping is reduced to half-filling, we find that 
quasiparticle (QP) dispersion becomes essentially unrenormalized, but in sharp contrast the QP 
spectral weight renormalizes to nearly zero. This opposite doping evolution of the renormalization 
factors for QP dispersion and spectral weight conspires in such a way that the specific heat remains 
Fermi liquid like at all dopings in accord with experiments. 

PACS numbers: 71.10.Hf,71.18.+y,74.40.-n,74.72.Kf 



I. INTRODUCTION 

Understanding how 'non-Fermi-liquid' behavior arises 
near the half-filled insulating state is one of the key ques- 
tions for unraveling the physics of not only the cuprates 
but that of correlated electron systems more generally. 
Here we show how some aspects of the electronic spec- 
trum which are difficult to understand in a conventional 
Fermi liquid theory can be explained naturally in a model 
of an antiferromagnetic Fermi liquid. Specifically, it has 
been reported that the renormalization of the electronic 
dispersion decreases with underdoping as half-filling is 
approached,— i.e. the quasiparticles (QPs) seem to 'un- 
dress' in the underdoped regime. In sharp contrast, the 
spectral weight of the QPs fades away on approaching the 
insulator and renormalizes to zero at half-filling^''^, indi- 
cating that these QPs are very fragile or 'gossamer'-like^. 
Within Fermi liquid theory the QP dispersion and spec- 
tral weight should be renormalized by the same factor, 
unless the self-energy has a strong momentum depen- 
dence. However, fluctuations in the cuprates seem to be 
relatively isotropic^, which would imply then that the 
renormalization factor Zd for dispersion is roughly equal 
to the renormalization factor for the spectral weight. 
This is clearly violated near the Mott insulating limit, 
where Zd ^ 1 while — >• 0. Added to these puzzling 
findings is the fact that the electronic specific heat con- 
tinues to behave more or less in a Fermi-liquid manner 
over the entire doping range from the overdoped metal to 
the insulator— i2r—. These results clearly demonstrate that 
a non-Fermi-liquid or 'strange metal' superconductor 
emerges from the Fermi-liquid background as doping is 
reduced. Notably, there is considerable controversy over 
what constitutes non-Fermi-liquid behavior and many 
proposals have been made to understand its possible ori- 
gin, ranging from preformed d— wave pairs^S to fluctu- 
ating spin-density wavesii, from the Hubbardi^ii^ and 
t — Ji"^ models to the anti-de Sitter/conformal field the- 
ory (AdS-CFT) correspondencei^. 

Insight into the origin of non-Fermi-liquid behav- 



ior comes from recent angle-resolved photoemission 
spectroscopy (ARPES)^^, Hall eSeci^, and quantum 
oscillation^^"™ experiments, which reveal the presence of 
FS reconstruction as a robust feature of both electron and 
hole doped cuprates, suggesting that the ground state in- 
volves some superlattice order. In this article, we investi- 
gate a model of the underdoped cuprates with a density 
wave ordered ground state and find that when the Fermi 
surface (FS) breaks into pockets, the self-energy devel- 
ops a strong momentum dependence. Analytic forms for 
various renormalization factors are presented to delineate 
how the non- Fermi-liquid physics arising from FS recon- 
structions can be understood at a quantitative level. Our 
study thus provides a tangible model for reconciling the 
seemingly contradictory doping evolutions of the QP dis- 
persion, the QP spectral weight and the specific heat in 
the cuprates. 

The possibility of FS reconstruction in the cuprates 
has been proposed many times since it was found that 
the Hall density scales with doping^! - i.e., it is pro- 
portional to the small pocket area rather than the large 
FS. We model this FS reconstruction by assuming a spin 
density wave (SDW) order which results in a nearly- 
antiferromagnetic-Fermi-liquid (NAFL) phase^^. In ad- 
dition to SDW order, we have analyzed other candi- 
dates for the competing order including charge, fiux, or 
d— density waves^'^ and find that the pseudogap symme- 
try which is the essential ingredient for the origin of non- 
Fermi-liquid physics is insensitive to the particular nature 
of the competing order state. We calculate the self-energy 
due to spin and charge fiuctuations within a QP-GW for- 
malism described in Appendix A. The model is known to 
properly describe the high-energy 'waterfall' features in 
ARPES^i^i2^ as well as the doping evolution of the opti- 
cal spectra and the finite QP lifetime. Within the QP- 
GW model the self-energy splits the spectrum into co- 
herent in-gap states and an incoherent residue of the un- 
doped upper and lower Hubbard bands (U/LHBs). With 
underdoping, the in-gap states develop an SDW gap and 
split into upper and lower magnetic bands (U/LMBs)^, 
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and the resulting 'four-band' features are consistent with 
quantum Monte Carlo (QMC) calculations^! [Appendix 
A]. The involvement of a quantum critical point (QCP) 
in the optimal doping regime where the SDW order dis- 
appears is strongly suggestive of an intermediate strength 
for correlations in the cuprates.— Accordingly, there have 
been several recent attempts to develop an intermedi- 
ate coupling model for the cuprates starting from either 
the weak-coupling limit, as in the present calculation, 
or the strong-coupling limit^^. Therefore, we compare 
our results with experiments as well as with the varia- 
tional cluster calculations of Paramekanti, Randeria and 
Trevedi (PRT)^^, which approach the problem from the 
RVB limit. 

This article is organized as follows. In sections IIA and 
IIB, we investigate the doping evolution of the spectral 
weight and that of dispersion renormalization, respec- 
tively. The specific heat results are presented in Section 
III. The discussion and conclusions are given in Sections 
IV and V respectively. A summary of the underlying QP- 
GW model is provided in Appendix A. Various renormal- 
ization factors invoked are summarized in Appendix B for 
the reader's convenience. Appendices C and D address 
details of analytic forms for the renormalization factors 
and those of our specific heat evaluation. 

II. NEARLY- ANTIFERROMAGNETIC-FERMI- 
LIQUID RENORMALIZTION 
FACTORS 

A. Momentum density and spectral weight 
renormalization 

We evaluate the exact values of the renormalization 
factors from FS discontinuities of spectral function mo- 
ments A/;(k) of various order /, 

/oo 
dojij^A{k,u)f{uj) (1) 
-oo 

These moments provide important information about the 
spectral weight distribution in energy and momentum 
space as a function of doping,^^ The spectral density 
A(k, uj) involves both the coherent QP and the incoherent 
part. Due to the Fermi function f{uj), the moments Mi 
display singularities at the Fermi momentum kp which 
are characteristic of coherent gapless quasiparticle exci- 
tations. We consider first the zeroth order moment which 
is simply the momentum density n(k)i^^— The interact- 
ing FS is determined from the jump in n(k) at kp of the 
quasiparticles. The magnitude of this jump defines the 
spectral weight renormalization Zi^ . The incoherent part 
in the spectral weight substantially modifies the shape of 
nik). 

Fig. [1] shows maps of n(k) throughout the first Bril- 
louin zone as a function of doping for Nd2-£!;Cea:Cu04 
(NCCO) and hole doped La2-a;Sr3;Cu04 (LSCO). In the 
present NAFL case the combination of self-energy and 
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FIG. 1: (Color online) Momentum density, n(k) as calculated 
using Eq. [1] is shown for various dopings for NCCO in (a-c) 
and for LSCO in (d-f). 



SDW coherence factors leads to characteristic structures 
in n{k) at all dopings including half-filling. At half- 
filling n(k) shows a maximum at the F— point and away 
from that it decreases gradually and smoothly, from in- 
side to outside the LDA-like FS [magenta solid line in 
Fig. [IJa) and (d)]. As we dope the system with elec- 
trons, the spectral weight increases at the F— point and 
in addition, (tt, 0) and its equivalent points largely gain 
spectral weight due to the development of electron pock- 
ets in NCCO [Fig. [Ijb)]. With further increase of dop- 
ing, the FS undergoes two topological transitions'^^''^^ as 
also refiected in the momentum density calculations here. 
The first topological transition in n(k) occurs when the 
LMB approaches the Fermi level (Ep) and forms hole- 
like pockets at (±7r/2, ±7r/2). For x = 0.15, the hole 
pockets are fully formed and they as well as the electron 
pockets increase in size with further doping. At a; = 0.18 
the electron and hole pockets merge at the hot-spot, the 
SDW gap collapses, and the full metal-like n(k) appears 
[second topological transition]. For hole doping, the FS 
topological transition is complimentary to the electron 
doped one and here the hole pocket appears first as shown 
in Fig. [Ije) and above the QCP, the electron-like full FS 
appears in Fig. [Ijf). 

A more quantitative account of the effect of self-energy 
corrections on the residual coherent QP spectral weight 
is provided in Fig. [2ja), which shows n(k) along high- 
symmetry lines for NCCO as well as LSCO. Some im- 
portant effects of correlations on the insulating state 
should be noted here. At a; = 0, n[k = F] « 0.9 and 
n[k = (tt,?:)] w 0.1, implying that the self-energy redis- 
tributes the spectral weight from the filled states to the 
unfilled regions even in the insulating phase. PRT [with 
a different parameter set t = 300 meV, t' = —t/A and 
U = I2t] find a similar result n\k = F] « 0.852^. At 
half-filling n(k) is a smoothly varying function through- 
out the BZ, due to the absence of gapless quasiparti- 
cle for both NCCO (green line) and LSCO (not shown). 
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As we increase electron doping, the spectral weight at F 
[(tt, tt)] gradually increases [decreases] whereas the spec- 
tral weight increases rapidly at (tt, 0), due to the devel- 
opment of electron pockets, while discontinuities in n(k) 
arise at the Fermi surface. In the underdoped region, 
n(k) shows additional singularities along F (tt, 0)) and 
(7r/2, 7r/2) — > (tt, tt) due to the presence of shadow bands 
as marked by gold arrows. Since the shadow bands are 
usually weak in the cuprates^"— , experimental data are 
available predominantly along the arcs - that is, along 
the antinodal direction for NCCO and nodal direction for 
LSCO. These are compared with our theory in Fig.[21[b), 
effects of the ARPES matrix element notwithstanding!^ 
We define a coherent spectral weight renormalization 
factor from the discontinuities in n{k): 



Anikp), 



(2) 



plotted as a function of doping for NCCO in Fig. [^Jb)^*. 
For the main band in NCCO along the antinodal direc- 
tion, the UMB crosses the Fermi level at all finite dopings 
(electron pocket) and Z^^ decreases smoothly with un- 
derdoping but vanishes discontinuously at a; = 0. Shown 
also in Fig.[21[b) is the conventional paramagnetic Green's 
function renormalization factor 
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It can be seen that the doping dependence in is very 
weak and strikingly opposite to Z^. Thus the strong 
doping dependence of Z^ is governed by the SDW gap 
collapse. In Appendix C we derive an approximate ana- 
lytical form for the SDW corrected renormalization factor 
(at w = 0) in terms of the SDW coherence factors 
[see Eq. [C6] 



^SDW _ fiu 
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Note however that kp is the true Fermi momentum in 
the SDW state, not the LDA one. We see in Fig djb) 
that this simple analytic form captures the essential dop- 
ing dependence of the full As doping increases, the 
weight of the shadow bands (above the magnetic Bril- 
louin zone) decreases [open symbols in Fig. H^b)], with 
the corresponding spectral weight shifting to the main 
bands. Along the nodal direction, remains zero up to 
optimal doping, and then shows a jump to a slowly in- 
creasing value as the LMB crosses the Fermi level. These 
results are in excellent agreement with experiment^. 

A similar doping dependence of Z,^ is observed in 
LSCO along the nodal direction, including the jump at 
X — and the spectral weight transfer from the shadow 
band to the main band^inFig.IlJb). This is in qualita- 
tive agreement with the calculation of PRT^^ (triangles 
in Fig.[2Ib)), and with ARPES measurements on LSCO^ 
(open circles) above the optimal doping region. However, 



in the very underdoped region, the ARPES data seem to 
extrapolate smoothly to zero at half- filling (dashed line) , 
which may be related to nanoscale phase separations be- 
lieved to be significant in LSCO. 

Note that the analytic renormalization factor Z^^^ is 
defined in the SDW phase, but is treated as a renormal- 
ization of the paramagnetic phase dispersion. This subtle 
point is an attempt to treat the pseudogap physics, and 
is discussed in more detail in Section IV below. 



B. FERMI VELOCITY AND DISPERSION 
RENORMALIZATION 

We turn next to the first order spectral moment Afi(k). 
One can measure the dispersion renormalization from the 
size of the slope discontinuity, which can be written as 
(Eq.Ell 



A{dMi{k)/dk)^=Z^VF, 



(5) 



where vp is the Fermi velocity^. Knowing Z^ from 
Figs. I^b) and (c), we can extract vp as a, function 
of doping, as seen in Fig. ^a.). The results are ob- 
tained for both LSCO and NCCO and compared with 
ARPES data on LSCO^ and with PRT's variational 
cluster calculations^'^ (hole doping) as well as with LDA 
and mean-field theory for LSCO. Mean- field results for 
NCCO have an equivalent doping dependence and thus 
are not shown here. Notably, in spite of the substan- 
tial decrease of Zi^, the QP velocity vp does not dimin- 
ish on entering into the pseudogap region. Although in 
the SDW mean field case vp decreases smoothly with 
underdoping as the gap grows, when a self-energy is in- 
troduced, the reduction of the coherent spectral weight 
(Zi^) with underdoping compensates for this, leading to a 
net enhancement of vp. The results are consistent with 
PRT and ARPES. Similar results are also obtained in 
a self-consistent Born approximation with antiferromag- 
netic pseudogap in the t — J model'*^ . 

We define a dispersion renormalization factor in terms 
of the conventional band velocity with respect to its bare 
(LDA) counterpart (vp) as 



Zd = vp/vp, 



(6) 



plotted in Figs. [3l^a). In the paramagnetic case, disper- 
sion renormalization is defined as 



_ .7O _ 



I + dJ:'/v"pdk) 



(7) 



This implies that in a conventional picture the devia- 
tion in behavior of Z^ from spectral weight renormal- 
ization Z^ comes solely from the A:— dependence of the 
self-energy. Since this self energy is approximately k- 
independent, one finds Z^ w Z^. In sharp contrast, the 
calculated Zd [Fig.[3Ib)] shows a strikingly opposite dop- 
ing dependence to Z^j. This can be understood to be 
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FIG. 2: (color online) (a) n(k) for NCCO and LSCO at several representative dopings with dashed lines showing the discontin- 
uous jump ai kp (highlighted in inset). Gold arrows mark features from the shadow bands which cross Ep- (b) Z^i at Ep are 
shown along the antinodal (red) and nodal (blue) direction for NCCO. Filled (open) symbols give the main (shadow) bands, 
compared with their corresponding analytical approximation Z^^^ plotted by solid (dashed) line of same color. ARPES result 
(green) is extracted from Ref. [l^. (c) Same as (b) but for LSCO along the nodal direction. These results are compared with 
PRT's calculations for hole doping2^ and ARPES results^ for LSCO along the nodal direction. All the experimental data and 
PRT data in (b) and (c) are normalized to highlight their doping evolution. Brown dashed line shows that if there is nanoscale 
phase separation in LSCO then Zu, would scale linearly with doping in the extreme underdoped region. 
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FIG. 3: (color online) (a) Fermi velocity vp along the antinodal direction in NCCO and nodal direction in LSCO is compared 
with ARPES data on LSCO?''"' and with PRT's calculation^-. The blue dashed (solid) line gives the corresponding LDA (MFT) 
results for LSCO. (b) Dispersion renormalization Zd for NCCO along antinodal direction, compared with an approximate 
analytical formula for the dispersion renormalization Z^^^ (solid line), (c) Same as (b) but for LSCO along the nodal 
direction. The results are compared with experimental dataii. 



the result of an SDW gap, which introduces a new fc- 
dependence in the dispersion renormalization as given 
by (Eq.EH): 

= ZIZIS^ = Z^Jl + -^-) . (8) 

Figs. EJb) and (c) compare Zd and Z^^^ for NCCO 
and LSCO respectively. Thus the doping dependence of 
Zd implies that as we go towards the Mott insulator, the 
dispersion tends towards the LDA-bands, consistent with 
LSCO results (blue open circles )i. The opposite doping 
dependences of Zd and Z^ can be readily understood by 
comparing the corresponding analytical formulas Z^^^ 



(Eq.ig) and Z^^^ {Eq.\^. Zf^^' {Z^^^) varies with 
SDW gap as A (1/A), increasing (decreasing) with un- 
derdoping. This in turn moves the in-gap states further 
away from the Fermi level (hence decreasing Z^^)^ and 
thus shifting the band towards the LDA values (increas- 
ing Zd)^. 



III. SPECIFIC HEAT 

A striking application of the renormalization effects 
can be seen in the doping evolution of the specific heat. 
A general expression for the specific heat cy in the SDW 



5 



4 



1 

(a) NCCO 


1 1 




~- 












—>''' • • LDA - 








. QP-GVr 


1 


• ■ Expt. 

1 1 



10 



5 



0, 



1 1 
(b) LSCO 


1 1 1 




. _ n 




— ° s Expt. - 


^ 1 J 
P X / ° 


^ X 


jff ^ J 
//^ 

^ 1 1 


— QP-GW 

1 1 1 



X 



FIG. 4: (color online) (a) Specific heat coefficient 7(2;) for 
different theoretical calculations (various lines), and the ex- 
perimental results on NCCO (red squares)* and PLCCO 
(red circles)?^. (b) Same as (a) but for LSCO (red squares"' 
and triangles^). Results are compared with ris [filled green 
squares'' and open green squares^] and x (blue)*''), all nor- 
malized at the VHS. The red dashed line shows that in un- 
derdoped LSCO 7 scales linearly with doping. Theoretical 
7 has been scaled by a factor of 1.1, consistent with a weak 
electron-phonon renormalization. 



for X < 0.11, reflecting the quasi-two dimensionality of 
the electron pocket (constant DOS) in cuprates. The step 
at X ~ 0.11 signals the appearance of the hole pocket. 
Similarly in LSCO the appearance of the electron pocket 
near (tt, 0) at x ~ 0.17 within our model greatly enhances 
the VHS features. Finally, adding self-energy corrections 
(in the form of preserves the general shape of the 

SDW 7, -while shifting its magnitude back to-wards the 
LDA values. It is interesting to relate the doping de- 
pendence of 7(0) and in Fig. [2][b). Note that Z^^ 
is evaluated at a particular Fermi momentum whereas 
7(0) is computed after summing over Zi^ at all kp. As 
Zi^ exhibits complimentary doping dependencies for main 
and shadow bands, so the total remains fairly constant. 
These results are in striking contrast to the strong cou- 
pling limit where 7 should diverge with the effective mass 
asx^ 

The agreement with experiment is quite good in LSCO 
for X > 0.10, including the strong VHS feature. Inter- 
estingly, the superfluid density and the paramagnetic 
susceptibility data, which are proportional to the total 
FS areas, show a similar doping dependence to 7. How- 
ever, below X = 0.10 7 as X 0, an effect not 
captured by our calculations. Such an effect could be 
due to a Coulomb gapi and/or nanoscale phase separa- 
tion. The linear dashed line in Fig. 3(b) illustrates the 
corrected form expected in the latter case. Note that 
nanoscale phase separation would produce the dashed 
line seen in Fig. ^c) as well as explain the anomalous 
doping dependence of the chemical potential."*^ Hence 
the enhanced gossamer features seen in LSCO may be 
related to nanoscale phase separation. 



state including the self-energy correction is given in Ap- 
pendix D. We find that the Sommerfcld coefficient 7 can 
be well described in terms of the density of state (DOS) 
(7V(0)) in the SDW state as^, 

l = cv{T)/T^^N{0)/Zl (9) 

where A'o(O) is the mean- field DOS with SDW gap but 
without self energy corrections. 

Fig. m compares experimental values of 7 as a function 
of doping in LSCO^ and electron doped NCCO® and 
Pri_a;LaCea;Cu04 (PLCCO)S with several calculations 
including bare LDA, MFT results with SDW gap, and 
with self-energy correction for both NCCO and LSCO. 
The striking differences between NCCO and LSCO away 
from half-filling are due to the presence of the van-Hove 
singularity (VHS) near Ep in the latter case. N{0) 
in LDA thus decreases (from a finite value at x = 0) 
with increasing electron doping, whereas for LSCO N{0) 
has a peak at the doping corresponding to the VHS 
xvHS ^ 0.20^. SDW order opens a gap, reducing A^(0), 
and introduces steps associated with the collapse of the 
SDW gap. Thus, for electron doping iV(0) is nearly fiat 



IV. DISCUSSION 

The present calculations are most appropriate for elec- 
tron doping, where only the (tt, tt) commensurate SDW 
order is observed, and the model is in very good agree- 
ment with experiment. Remarkably, the same model 
when applied to hole-doped cuprates describes many as- 
pects of the two-gap scenari o^^'^-^ , despite the fact that 
it does not capture the incommensurate magnetization. 
The non-Fermi-liquid behavior presented here is not sen- 
sitive to the specifics of the competing order, but only 
to the resulting superlattice g-vector. We have ana- 
lyzed other candidates for the competing order includ- 
ing charge, flux, or d— density waves—, and find that the 
results are insensitive to the nature of the competing or- 
der state. Thus, Eqs. 4 and 8 continue to hold for any 
Q — (tt, tt) order, as long as the appropriate gap A is 
used. 

In the presence of long-range magnetic order, the 
Green's function develops a second pole, and should 
properly be treated as a tensor. However, the pseudogap 
phase is more likely to be associated with only short range 
order, in which case the second pole does not cross the 
real axis, the Green's function can be treated as a scalar. 



6 



and the information about incipient gap formation is en- 
coded in the self-energy^ The analytic approximations of 
Eqs. 4 and 8 mimic this effect by treating the renormal- 
ization factors as acting on the paramagnetic dispersion 
^k, with information on proximity to long-range magnetic 
order encoded into Z^^^ and Z^^^ . 

This approach is similar to the phenomenological 
model introduced by Yang et al.— to describe RVB 
physics. Indeed, their phenomenological self energy is 
quite similar to the form expected for a NAFL, except for 
the treatment of scattering at the magnetic zone bound- 
ary, which splits the pockets expected for AF order into 
two half-pockets. One puzzle is that the effect of strong 
scattering at a superlattice zone boundary is already 
well understood in standard Lifshitz-Kosevich theory^, 
where it leads to a quantum switching between arcs of 
Fermi surface. This produces harmonic mixing frequen- 
cies in a quantum oscillation (QO) spectrum, and it is 
hard to see how it could evolve into the short-circuiting 
effect of Yang, et al. At any rate, the QOs observed in 
NCCO are more consistent with the present modelJ^ 



V. CONCLUSION 

In conclusion, we have shown that a number of salient 
features of the non- Fermi-liquid state of the underdoped 
cuprates can be understood within the framework of a 
competing density wave order, which breaks the particle- 
hole symmetry, and drives reconstruction of the FS. We 
provide a transparent and analytic basis for describing 
how the non-Fermi-liquid effects play out in renormal- 
izing spectral weight via Zi^ and electronic dispersion 
via Zd, and how they conspire to yield a specific heat 
in the cuprates which is essentially conventional in na- 
ture at all dopings. Our framework would provide a 
straightforward basis for understanding how the broken- 
symmetry order leads to 'non-Fermi-liquid' effects, not 
only in the cuprates, but also in heavy-fermions^, Fe- 
based superconductors^^ and other strongly correlated 
materials. 
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Appendix A: Details of quasiparticle-GW (QP-GW) 
Model 



In the QP-GW formalism the bare dispersion is taken 
as the LDA dispersion (^k = Ck — -^'-f): modeled via a 
tight-binding (TB) fitS^H. We calculate the self-energy 
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FIG. 5: (color online) (a) Average renormalization factor Z 
decreases linearly with doping, very much like in Figs.[2l^b) 
and (c). (b) Our computed doping dependence of selfconsis- 
tent values of U and ZU is compared with earlier mean field 
results^i^. 
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FIG. 6: (color online) (a) Spectral intensity as a function of 
oj along the high-symmetry lines for several dopings (at tem- 
perature r = 0). Blue to red color map gives the minimum 
to maximum intensity. The yellow dashed line gives the un- 
derlying LDA dispersion where the gold lines represent the 
renormalized magnetic bands (So-dressed), (b) The QP-GW 
DOS (blue lines) is compared with Eo-dressed DOS (red line) 
calculated at T = 0. The green lines show the DOS without 
the vertex correction. 



in a GW-like formalism using a simplified (one param- 
eter) scheme where the input and final self-energies are 
self-consistent in the coherent part only (QP-GW model). 

In the following, we use a 'tilde' over a quantity to 
symbolize that it is a 2 x 2 matrix. The self-energy in 
the underdoped region is written in canonical form using 
the Nambu formalism as 

i;*(k, cr, iiLj„) = i;(k, CT, iwn) + (/)(k, cr, iw„)fi. (Al) 

The dressed Green's function is then 

G~-^(k, CT, iw„) = iujnl - Hlda - S*(k, cr, ia;„) (A2) 

where H^oA is the bare Hamiltonian defined in the 
magnetic zone as diag[^k^, ^k+q], and the fi are Pauli 
matrices. At each step, we adjust the chemical potential 
Ep to fix the doping. Real frequency Green's functions 
are extracted from the Matsubara results by analytic 
continuation iujn — ^ uj -\- iS. Here, S](k, cr, «aj„) is a 
2x2 matrix in the SDW stated, whose diagonal part 
renormalizes Hlda, while the off-diagonal term gives 
a (small) anomalous frequency dependence to the gap 
function (/)(k, cr, ia;„) = aA. A = US is the magnetic 
gap parameter and S is the magnetization at the 
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commensurate vector Q = (7r,7r) for Hubbard C/, which 
is calculated using a mean-field approximation^^. The 
doping dependence of the on-site Hubbard U is obtained 
due to charge screening from U — {V{q)/{1 + V{q)x{q)), 
where V{q) is the long-range Coulomb interaction^, 
and xil) is the charge susceptibility in the gapped 
state defined below. The obtained values of screened U 
are given in Ref. [H and plotted in Fig. [HJb). 

We calculate the self energy due to the spin as 
well as the charge response within a GW framework as 

S(k,a,za.„) = |L/^Z^ 

G(k + q, a', iw„ + Wp)r(k, q, iu„, Wp)Im[xRpA(q, Up)], 

(A3) 

where the prime over the momentum summation indi- 
cates that the summation is restricted within the mag- 
netic Brillouin zone. The dressed susceptibility in the 
above equation is given in terms of the 2x2 RPA sus- 
ceptibility asr^ 



susceptibilities become 



xrpa (q>«^n) 



(A4) 



Here the superscript iaa) refers to the combined charge 
plus longitudinal spin susceptibility tensor, whereas (era') 
(with a — — fj) gives the transverse susceptibility tensor, 
and the Xq'^^'^'^ (q, jcj„) are bare susceptibilities. 



A self-consistent 'dressed' GW calculation would 
include S on the right-hand side of Eq. IA3I in both 
G and xo- This generally leads to problems unless a 
vertex correction F is included. Improved results are 
often found by using a 'bare' Go and xo ~ bare in the 
sense of not including E. This latter GoWq— scheme 
also fails in the present calculation by producing too 
large a renormalization of the band dispersion. This is 
because the imaginary part of the bare susceptibility 
has the form Xo ~ ^(^^ ~ ~ '?fc])j so that near 

the Fermi surface, Xo should scale in frequency with 
the dressed quasiparticle dispersion. Since Go Wo uses 
the bare dispersions, peaks in Xoi which control the 
renormalization, lie at too high an energy. 

We therefore introduce a modified, or 'quasiparticle' 
(QP) GW approximation^ for the right-hand side (RHS) 
of Eq. IA3[ as follows. In Eq. IA3[ we dress both Go and 
Xo with an 'input' self-energy chosen as 



E(za;„) = (1 - 1/Z)iw„i. 



(A5) 



Thus, the input E contains a single parameter Z, which 
gives an overall renormalization of the 'input' dispersions 
(RHS of Eq. IA3|) . With this approximation, the bare 



~aa /(T(T 

Xo 



(q, = -z'yy sz/^^[^^5^l_I^^<l. 



k v,v' 



(A6) 



The prime over the summation has the same meaning as 
in Eq. IA3I In Eq. IA6[ the v summation is over the two 
magnetic bands UMB [v = +) and LMB (zy = -): 



(A7) 



with ^± = (^k ± a+Q)/2 and Eoi. = ^ {^^ ^ + 
(Ref. 0), f{E) = 1/(1 + exp{E/kBT)) is the Fermi 
function at temperature T, and fc^ is the Boltzmann con- 
stant. The coherence factors S'^'^J,"'^ give the amplitude 
of the scattering of the quasiparticles with the charge 
and magnon modes of the system respectively with com- 
ponents 

^^.^"''(ll) = (ak«k+q ± ^^^^'/3k/3k+q)^ 



'(7(7 1 (J(7 



(12) 



"i^(ak/3k ± i^i^'ak+q/3k+q)- (A8) 



Here 



(A9) 



respectively are the weights associated with the U/LMBs. 
The other coherence factors in Eq. IA8I can be derived 
using the translational symmetry with respect to q. 

A limitation of the present scheme is also apparent 
from Eq. IA3I This self-energy is a good approximation 
for the coherent, dressed bands, but does not extend to 
the incoherent part of the spectrum, thereby underesti- 
mating the incoherent spectral weight. We have empir- 
ically found that this can be partly remedied by incor- 
porating a vertex function F. Consistent with the QP 
approximation, the vertex correction to the self-energy is 
obtained through Ward's identity as 

F(k,q,iw„,wp) = 1 - (aE79t^)„=^„ = 1/Z.(A10) 

Even though we have chosen perhaps the simplest form of 
the vertex correction, it has a large impact on the spectral 
weight transfer. F = 1/Z eventually reduces the renor- 
malization of the bare susceptibility in Eg. IA3l and hence 
the spectral weight is spread out more towards higher 
energies, enhancing the incoherent spectral weight. An 
illustration of the importance of this vertex correction is 
given in Fig. [HKb). 

Next we discuss our self-consistent scheme, and ex- 
plain how the parameter Z is chosen. We choose Z to 
match the average renormalization in the low-energy (co- 
herent) part of the spectrum. Specifically, if E' is the 
real part of the diagonal self energy, then we adjust Z 
self-consistently until it satisfies Z = {1 — dJ^' / dui)^^^ , 
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FIG. 7: (color online) Real and imaginary part of the self- 
energy as expanded in the tight-binding form of Eq. All. 



where Wo is an average quasiparticle excitation energy, 
which is related to the poles in G. This gives a good self- 
consistent result for the coherent spectral weight in the 
low-energy region (see Fig.[6][a)), whereas the incoherent 
parts in the higher energy regions are not self-consistent. 
Thus our scheme is in the spirit of Landau's quasiparti- 
cles, except that Landau assumed that all of the spectral 
weight goes into the QP band, while we have only a frac- 
tion Z . 

When this is done, we find that U is effectively further 
renormalized by the doping dependent Z [Fig. [Sja)], so 
that the product Uxo is approximately independent of Z . 
The resulting ZU closely resembles our earlier mean-field 
calculations as shown in Fig. [Sljb). 

Lastly, we find that the momentum dependence of the 
fluctuation self-energy S is relatively weak^'^. To empha- 
size this point, we expanded the momentum dependence 
of the self energy in a form similar to the tight-binding 
model: 

E(i:j,k) = T.(3[u) + (cx + Cy) 

+ T,2{uj)c.j,Cy + S3(w) {C2x + C2y) ,(All) 

where Ca(x/y) = cos (Q;fc(^/j,)o). We calculate the 
self-energy at four high symmetry points k = 
(0,0), (tTjO), (7r,7r) and (7r/2,7r/2) to obtain the above 
coefficients as shown in Fig. [T] Clearly, only the 
fc— independent part (red line) has a strong contribution. 
Therefore, we have simplified the self-energy calculation 
by approximating it with a fc-independent average value 
taken as the value at fc = (7r/2,7r/2). Since we are ne- 
glecting the /c— dependence of the self energy. En — E22 
and E12 = E21. 

FigurelSJa) shows that in the QP-GW scheme the spec- 
tral weight splits up into four band-like features. The two 
bands closer to the Fermi level are the UMB and LMB, 
associated with the development of the spin-density wave 
(SDW). The residual incoherent spectral dispersions at 
higher energies are the UHB and LHB. Similar four band 
features are found in the variational cluster calculations 
shown in Fig.[51Jb)2I. The corresponding density of states 
(DOS) shows four peaks associated with these four bands. 
Figs. [Hc)-(f). These four band features separated by 
SDW gap or 'waterfall' effects show semi-quantitative 
agreement with QMC result o^'^^ . With doping, the two 
magnetic bands merge in a QCP near optimal doping. 



(a) QP-GW (This work) (b) Variational Calculation 



x=0.0 



ABPES, T - O.IO; SDW fit — 




r (71,0) (7t,7t) r r (ti.O) (ti.ti) r 

k[1/a] 



QMC 



Present calculations 




-2 2 
ffl (eV) 



2 2 
ffl (eV) 



FIG. 8: (color online) (a) The spectral intensity of NCCO at 
a; = 0, plotted (in logarithmic scale) along the high-symmetry 
lines, is compared with variational calculations^^ in (b). (c)- 
(f) The computed DOS at various dopings are compared 
with the corresponding QMC results (blue lines) for x = 0.0 
[Ref.im and for x = 0.05 to x = 0.20 [Ref.iej. 



while the two Hubbard bands occur at the top and bot- 
tom of the LDA bands, and hence the associated Hub- 
bard band splitting is comparable to the LDA bandwidth 
at all dopings. 



Appendix B: Summary of renormalization factors 

We summarize the various renormalization factors that 
arise in this study for convenience reference. Equation 
numbers in square brackets indicate where a specific fac- 
tor is discussed in the text. 

Z^ = Anikp) [Eq.[2] 
— Total spectral weight renormalization. 

Z° = (1 - 9S79c.):i„ [Eq.E 
= Self energy contribution to spectral weight 
renormalization. 
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tSDW 



ySDW 



2 



1± 1 



2A 




-1/2- 



[Eq.i] 

Analytical formula for SDW spectral weight 
renormalization. 



Zd = vp/vp 

= Velocity renormalization. 



[Eq.[6] 



= z°^zl^zl{i + di:'/v%dk) [Eq.[7] 

— Conventional dispersion renormalization. 



(1 + AVa.a.+g) [Eq.lH] 
Analytical formula for SDW dispersion 
renormalization. 



where £",, 



Z'^.E^] 



Therefore, the coherent part is 



governed by the SDW coherence factors ak and /3k for 
the filled state with a renormalization by Z^, ■ The zeroth 
and first order moment of the spectral weight (Eq. [1]) 
then can be approximated respectively as^, 



n(k) = 



Mi(k) 



Aii(k, uj)duj 

Aii(k, Lj)ujduj 

p)[alv+6{-E+)+Plv-t 



(C4) 



-e: 



(C5) 



In the second equation above, we have expanded the 
renormalized band near the Fermi level as E^ w (k — 
ki?), where vp is the corresponding Fermi velocity. The 
singularity in n(k) becomes 



7 = cv/T 

= Specific heat coefficient. 



[Eq.[9] 



Z^ = An(k^^) = Zlal^ Z^ 



for UMB 

for LMB. (C6) 



Appendix C: Analytical form of various 
renormalization factors 



Near the Fermi level, the dressed Green's functions 
from Eq. IA2l can be approximated as. 



Gii(k,w) 



k+Q 



Z" 



ySDW 
^d 



(k)6 



(CI) 



where = Z°£,i^, A = Z°A, and ^k is the bare (LDA) 
dispersion calculated at the same doping. In Eq. ICll 
Z° = (1 -5S;i(w)/9a;)-i (Eq.© is the part of the spec- 
tral weight renormalization due to the self-energy, where 
we have neglected the contribution of E12. Even though 
we assume a fc— independent self-energy and thus ^k = 1 
from Eq. ICll the SDW gap introduces a fc— dependent 
dispersion renormalization (defined here at a; = 0) given 

by 



zSDW 



(k) = ZM 1 



A2 



(C2) 



Note that ^k, ^k-i-q, and A are aU bare values in Ea. lC2l 
as the same renormalization factor Z^ gets cancelled out 
in the last term. Furthermore, we can split the spectral 
function ^(k, w) = ImG(k, a;)/7r into a coherent part (at 
Fermi level) and an incoherent part where the former can 
be represented by a delta function as 



An(k,a.) = Z"^{alS{uj~E+)+l3lS{Lo^E^))^ 

~\~ -^incoh 

(k,LL; > 0) 



(C3) 



Note that in the present case, the weight for UMB (LMB) 
appears along the antinodal (nodal) direction when the 
band crosses the Fermi level. Inserting the form of a, /3 
from Eq. IA91 we get Eq. 21 Similarly, inserting this Z^ 
in Eq. IC5[ we can measure the singular jump in Afi(k) 
as 



A(dAfi(k)/dk) 



•^u ^P- 



(C7) 



Thus and vp acquire a ^-dependence through the 
SDW coherence factor. 



Appendix D: Specific Heat Calculation 

Following the derivation by Abrikosov et. al^, we 
calculate the entropy in the SDW state for a strongly 
correlated system at finite temperature (in the low tem- 
perature limit) as 



S{T) = 



ml 

2711 



dujLU — 



dm 

did 



k,cr 



Tr [in G]^^ (k, a, u) - In G^\k, a, u) 



(Dl) 



Here the 2x2 retarded and advanced (dressed) Green's 
functions Gr/Ga depend on the temperature through 
the SDW order parameter only, which has a very weak de- 
pendence in the low temperature region. We can rewrite 
Eg. lDll in terms of a dimensionless parameter y = /3a; fol- 
lowing Ref. [g^ . and then taking the temperature deriva- 
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tive we get the expression for the specific heat as 
cv{T)^-^ dyy^ sec\i\y/2) 



An 



Im ( G_R,(k, a,uJ)-^GJ^{^^,a,i^J] 



(D2) 



where An is given in Eq. lC3l and A22 = ^ii(k — >• k+Q), 
i.e., A22 is similar to An, only the weights (a^, /3^) are 
interchanged. We have used Eq. ID3I in the calculation 
of Fig. 5. The constraint + /^k = 1 removes these 
SDW coherence factors from the equation. Neglecting 
the off-diagonal term proportional to the small quantity 
dT^i2/duj, the specific heat expression becomes 



cy(T) behaves linearly with T in the low temperature 
region, while the slope (7) undergoes an abrupt change 
with a kink in the waterfall region. Since the high energy 
kink energies for cuprates are around 0.3 to 0.6 eV, the 
kink should appear in cv only at a very high temperature, 
Tk ~ 10'^ if. Thus, in our calculation of 7 (Fig. 5), we 
have used the full expression for cy above but evaluated 
it in the co ^ limit. 

In the T — {uj ^ 0) hmit, Eq. ID2l can be simplified 
as sech(y) = 5(y — /3a;). In this limit, the imaginary 
part of the self-energy is zero and thus the last quantity 
in Eq. ID2l is calculated by noting that the w— derivative 
of the real part of the self-energy is: diagonal term = 
1 — l/Z^, and off-diagonal term —{dY,[2/9i^)Lj=()- This 
simplifies Eq. ID2I as 



Cv 



1 



(An(k,0)+A22(k,0)) 



dui 



(Ai2(k,0)-|-A2l(k,0)) 



(D3) 



cv « T:n'klTY,[S{-Z'^E+) + S{-Z"^E^)] 



-n'klT [iV+(0) + iV-(0)] 



(D4) 



Here, 7V^(0) = J2k^i~^k) ^^"^ total density of states 
for the (U/L)MBs at the Fermi level in the SDW state, 
but without renormalization (mean- field theory). For 
both NCCO and LSCO, Eq. lD4l is an excellent approxi- 
mation to the exact expression of Eq. ID3I It is interesting 
to observe that in the final expression for cv (Eq. ID4p , 
the self energy correction enters only through the renor- 
malized DOS. Thus, in the QP-GW model, cy has the 
same form as in conventional Fermi-liquid theory. 
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